Precise comparison of theory and new 
experiment for the Casimir force leads to 
stronger constraints on thermal quantum 
effects and long-range interactions 
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j3 ■ Abstract 

We report an improved dynamic determination of the Casimir pressure p ex P* be- 
tween two plane plates obtained using a micromachined torsional oscillator. The 
main improvements in the current experiment are a significant suppression of the 
surface roughness of the Au layers deposited on the interacting surfaces, and a 
decrease by a factor of 1.7 (down to 0.6 nm) in the experimental error in the mea- 
surement of the absolute separation. A metrological analysis of all data for p ex Pt 
from 15 sets of measurements permitted us to determine both the random and 
systematic errors, and to find the total experimental error in p ex P* as a function 
of separation at the 95% confidence level. In contrast to all previous experiments 
on the Casimir effect, where a small relative error was achieved only at the short- 
est separation, our smallest experimental error (~ 0.5%) is achieved over a wide 
separation range. The theoretical Casimir pressures p theor in the experimental con- 
figuration were calculated by the use of four theoretical approaches suggested in the 
literature based on the Lifshitz formula at nonzero temperature. All corrections to 
the Casimir force due to grain structure of the overlying metal layers (including the 
variation of optical data and patch potentials), surface roughness (including non- 
multiplicative and diffraction-type effects), and nonlocal effects, were calculated or 
estimated. The maximum value of the roughness correction, achieved at the shortest 
separation of 160 nm, is equal to only 0.65% of the Casimir pressure. All theoretical 
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errors, including those introduced by the proximity force theorem, finite size of the 
plate area, and uncertainties in the experimental separations, were analyzed and 
metrologically combined to obtain the total theoretical error at the 95% confidence 
level. Finally, the confidence interval fbr (p theor — pexpt^ wag obtained as a function 
of separation. Our measurements are found to be consistent with two theoretical 
approaches utilizing the plasma model and the surface impedance over the entire 
measurement region from 160 nm to 750 nm. Two other approaches to the thermal 
Casimir force, utilizing the Drude model or a special prescription for the determi- 
nation of the zero-frequency contribution to the Lifshitz formula, are excluded on 
the basis of our measurements at the 99% and 95% confidence levels, respectively. 
Finally, constraints on Yukawa-type hypothetical interactions are strengthened by 
up to a factor of 20 in a wide interaction range. 
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1. Introduction 



Recent advances in the experimental and theoretical investigation of the Casimir 
effect suggest that this field of research is important for both fundamental 
and applied physics. The Casimir force, which acts between two uncharged 
plates in vacuum, is a purely quantum effect caused by the alteration of the 
zero-point oscillations of the electromagnètic field due to boundaries [1]. The 
modern stage in Casimir research originated in a number of measurements of 
the Casimir force between metals (see [2,3,4,5,6,7,8,9,10,11] and review [12]). 
Additional motivation for studying the Casimir force comes from promising 
applications of these results in both fundamental physics (for constraining hy- 
pothetical long-range interactions [11,13,14,15,16,17]), and in nanotechnology 
[18,19]. Recently the Casimir- Polder force has attracted attention [20,21] in 
connection with experiments on quantum reflection and Bose-Einstein con- 
densation. 

As often happens, experimental progress placed more stringent demands upon 
theoretical calculations of the Casimir force. At present, to achieve an agree- 
ment between theory and experiment at the level of 1% of the measured force, 
special care must be taken to account for all relevant factors, such as the fi- 
nite conductivity of the metallic boundaries, the effect of thermal corrections, 
surface roughness, grain structure of the covering layers, and patch potentials. 
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These are reviewed in [12], and more recently in [22] which is devoted to the 
application of these qüestions to Casimir force measurements between gold 
surfaces by means of an atòmic force microscope [6]. 

The bàsic theoretical formalism for the calculation of the van der Waals and 
Casimir forces between real materials is given by Lifshitz theory [23,24,25]. 
Beginning in 2000, several research groups applied this theory to the case 
of real metals at nonzero temperature and four different approaches were 
proposed. In the first approach [26] the zero-frequency contribution to the 
Lifshitz formula was determined by the use of the Drude dielectric func- 
tion. The second approach [27] modified the zero-frequency term of the Lif- 
shitz formula to make it the same for ideal and real metals. In the third 
approach [28,29] the Lifshitz formula was combined with the free electron 
plasma model neglecting relaxation. Finally, in [30] a fourth approach was pro- 
posed in which the Lifshitz formula at nonzero temperature is rederived using 
the Leontovich surface impedance boundary conditions rather than the bulk 
dielectric permittivity. Within this last approach [30,31], the zero-frequency 
term of the Lifshitz formula is determined by the form of the Leontovich 
impedance vàlid near the characteristic frequency. Detailed discussions of the 
advantages and disadvantages of each of proposed approaches can be found in 
[32,33,34,35,36,37,38,39,40,41,42]. These discussions are part of a broader con- 
text connected with the fundamental nature of the zero-point oscillations and 
their relation to bàsic physical principies such as the laws of thermodynamics. 
In particular, it has been shown [31,33] that the first and second theoretical 
approaches to the description of the thermal Casimir force between real met- 
als lead to a violation of the third law of thermodynamics (the Nernst heat 
theorem) in the case of perfect crystals. The attempt made in [37,38] to avoid 
this conclusion in the framework of the first approach required the introduc- 
tion of defects and impurities to the crystal, which cannot be considered as 
satisfactory since the perfect crystal forms the basis of the theory of electron- 
phonon interactions. It follows that any formalism applied to perfect crystals 
must satisfy the third law of thermodynamics, which the formulations of the 
first approach do not. 

We consider now a comparison of the different theoretical approaches with 
experiment. It should be noted that the third and fourth approaches predict 
small thermal corrections to the Casimir force at short separations of order 
100 nm, in qualitative agreement with the case of an ideal metal, whereas 
the first and second approaches lead to thermal corrections which are many 
times greater [32]. However, taking into account the small magnitude of the 
Casimir force, reliably distinguishing among the different predicted thermal 
corrections presents a serious challenge to experimental physics. In this regard 
it is important that the first approach (see [26] and also [34,35,36,37]) disagrees 
significantly [41,42] with the most famous modern measurement of the Casimir 
force between Au surfaces by means of a torsion pendulum [2] . A net deviation 
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from zero temperature Casimir force is predicted by the first approach to be 
about 25% at a separation a = 1000 nm. It was not, however, observed [2,41,42] 
despite the fact that the experimental uncertainty was less than 10%. 

The most precise and accurate measurement of the Casimir force to date was 
made between two dissimilar metals (Cu and Au) by means of a microelec- 
tromechanical torsional oscillator [10,11]. In that experiment the Casimir force 
between a plate and a sphere was measured statically, and the effective Casimir 
pressure between two parallel plates was determined dynamically, with an ab- 
solute error of rs 0.6 mPa at 95% confidence. The estimated theoretical uncer- 
tainty in [11] was approximately 1% of the Casimir pressure over the entire 
measurement range. This was, however, a purely theoretical error which did 
not take into account errors in the theoretical pressure induced by uncertain- 
ties in the experimental determination of the sphere-plate separations. (Recali 
that a micromachined oscillator was first used in [18] to demonstrate a mi- 
cromechanical device driven entirely by the Casimir force.) The experimental 
data in [11] were not precise enough to observe the traditional thermal cor- 
rections of the third and fourth approaches. However, the comparison of these 
data with theory demonstrates that the alternative thermal corrections to the 
Casimir pressure, as predicted in the first [26] and second [27] approaches, are 
excluded by experiment within a wide separation region. The good agreement 
between theory (based on the surface impedance approach) and experiment 
was used in [11] to obtain stronger constraints on hypothetical long-range 
interactions predicted by many extensions to the Standard Model of funda- 
mental interactions. 

In this paper we present the results of an improved experimental determination 
of the Casimir pressure reported in [11], and of a refined comparison of these 
new results with all four theoretical approaches. The following modifications 
in the experimental setup allowed us to obtain results exceeding those of [11] 
in both reliability and conclusiveness. The main improvement is associated 
with a considerable decrease of the surface roughness on both the plate and 
the sphere. In [11] the highest roughness peaks were as much as ~ 100 nm in 
height on both bodies. In the improved experiment the maximum heights of 
the roughness peaks are 11.06 nm on the sphere and 20.63 nm on the plate. 
An improvement in detection sensitivity, together with a reduction of the 
coupling between the micromachined oscillator and the environment, yielded 
measurements at smaller separations between the test bodies (160 nm instead 
of 260 nm). Another important improvement was the reduction of the error 
in determining the absolute separation between the sphere and the top of the 
metal-coated plate from lnm in [11] to 0.6 nm. A special experimental study 
was performed to ensure the linearity of the microelectromechanical torsional 
oscillator. In addition Au was used as a covering layer on both the plate and 
the sphere (instead of Cu and Au, respectively, as in [11]), which permits a 
more reliable theoretical interpretation of the experimental results. 
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The refined theory includes a complete analysis of the finite conductivity, 
surface roughness, and thermal corrections to the Casimir pressure. In con- 
trast to [11], the finite conductivity corrections are analyzed by taking into 
account the grain structure of the covering layers, possible sample-to-sample 
variations of the tabulated optical data, and the contribution of nonlocal ef- 
fects. The roughness analysis incorporates non-multiplicative, diffraction-type 
and correlation effects. It is shown that in the present experiment the maxi- 
mum value of the roughness correction, which is achieved at a separation of 
160 nm, is only 0.65% of the Casimir pressure. Thermal Casimir pressures are 
computed in the framework of all four theoretical approaches. 



A special metrological analysis has been performed separately for the experi- 
mental data and for the results of the theoretical computations. The random, 
systematic, and total experimental errors are found as functions of separation, 
taking into account the distribution laws for each error at the 95% confidence 
level. A careful comparison between experiment and theory has been carried 
out taking into account all experimental and theoretical errors. This compar- 
ison is based on rigorous metrological criteria which do not use the concept 
of the root-mean-square deviation between theory and experiment over the 
entire separation range, a procedure that had been criticized in literature [7]. 
As a result, the alternative thermal corrections to the Casimir pressure, as 
proposed in the first [26] and second [27] approaches, are excluded by exper- 
iment, whereas the third [28,29] and fourth [30,31] approaches are shown to 
be consistent with experiment. These results confirm those obtained earlier in 
[11]. In this paper we also improve the constraints on hypothetical long-range 
interactions obtained in [11]. By combining the results of this paper and [11], 
the previously known constraints are strengthened by a factor of up to 20 in 
a wide interaction region around 0.1 /ím. 



The present paper is organized as follows. In Section 2 the experimental con- 
figuration and sample characterization are briefly described. Section 3 contains 
the calibration procedure and the experimental results for the Casimir pres- 
sure. In Section 4 the metrological analysis of these results is performed and the 
experimental errors are found at 95% confidence. In Section 5 we calculate the 
Casimir pressure between Au layers in the framework of the four approaches 
described above taking into account different kinds of corrections. Section 6 is 
devoted to a metrological analysis of the theoretical uncertainties. In Section 7 
all four theoretical approaches to the thermal Casimir force are compared with 
the experimental data. The level of agreement between experiment and theory 
is used to obtain stronger constraints on hypothetical long-range interactions 
in Section 8. We conclude in Section 9 with a discussion of our results. 
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2. Improved experimental setup and sample characterization 

2.1. Details of the improved experiment 

The general scheme of the experimental setup used to determine the Casimir 
interaction between two bodies by means of a micromachined oscillator was 
described previously in [10,11]. One of the new features of the present ex- 
periment is the fact that here the Casimir attraction is measured between 
two layers of Au, specifically, a Au-coated sapphire sphere and a Au-coated 
polysilicon microelectromechanical torsional oscillator (MTO). To improve the 
force sensitivity, the vertical separation between the sphere and the plate was 
changed harmonically with time, leading to a measurement of the z-derivative 
of the Casimir force. As described in [11], this is equivalent to measuring the 
Casimir force per unit area, or the Casimir pressure, for a configuration of 
two parallel plates (see below). The MTO presents a low coupling with the 
environment, since it is less sensitive to modes that involve a displacement 
of the center of mass. The miniaturization also allows us to achieve a large 
quality factor Q ~ 8000. 

The experimental arrangement is shown schematically in Fig. 1. The MTO is 
made of a 3.5 /im thick, 500 x 500 /xm 2 heavily doped polysilicon plate sus- 
pended at two opposite points by serpentine springs. Two independently con- 
tacted polysilicon elèctrodes located under the plate are used to induce an 
oscillation on the plate at the required resonant frequency of the MTO. A 
80 /an wide ribbon at the edge of the plate was coated with 10 nm of Pt fol- 
lowed by 150 nm of Au. The Au layer constitutes one of the metals used in 
the measurement of the Casimir interaction. The separation between the two 
Au layers is given by 

z = z meas - D -b6. (1) 

In this expression z meas is the separation between the end of the cleaved fiber 
and the platform, as interferometrically measured with the absolute error 
Az meas = 0.2 nm. D = Di + D 2 , where D 2 represents the separation be- 
tween the platform and the top of the Au-coated plate when the interaction 
between the sphere and the plate is negligible (i.e., = 0), and D\ is the 
distance between the end of the bottom of the cleaved fiber and the external 
surface of the Au film on the Au-coated sphere (see Fig. 1). The lever arm 
b = (210 ±3) /iin is determined optically. The value of 9 is determined by mea- 
suring the difference in capacitance between the plate and the right and left 
elèctrodes C r ^\ = C r i g ht — Ci e ft. In all reported cases 9 < 10~ 5 rad. A field effect 
transistor (FET), the first stage in the amplifier circuit, was placed as close as 
possible to the oscillator to minimize the effect of parasitic capacitances. The 
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smallest angular deviation that can be detected in this configuration is given 
by 



Ad 



AC r j 2d g _ AV C T 2d g 



io- 9 v^ 



rad 



(2) 



where AC r j is the minimum detectable change in C T ^ and C Q ~ e ww'/d g is 
the capacitance between the oscillator and each elèctrode when no forces are 
present (e Q is the permittivity of free space). The dimensions of the elèctrodes 
are w = 250/iin, w' = 190/xm, and d g = 2/im is the separation between the 
plate and the elèctrodes. Ct ~ 20 pF ^> C Q is the parasitic capacitance of the 
measurement circuit (determined by the capacitance of the FET transistor 
plus parasitic capacitances to ground), AV ~ lOy^nV/Hz 1 / 2 is the input 
noise of the amplifier, v is the bandwidth of the measurement, and Vi (i = 
1,2) ~ 1 V is the DC potential used to correct for initial asymmetries in 
the circuit and to linearize its response. In obtaining Eq. (2) parallel plate 
approximations, and the fact that the maximum angular deviation 9 max <C 1, 
have been used. In the relevant (100 — 5000) Hz range, actual measurements 
of AC r j coincide within 15 % with those predicted by Eq. (2) . 

A sphere with a nominal radius R = 150 /xm was coated with a 10 nm layer 
of Ti followed by a 200 nm layer of Au. For both the plate of the MTO and 
the sphere, vacuum was not interrupted in between depositions. Deposition 
induced asymmetries in the sphere were found to be smaller than 10 nm, the 
resolution of the scanning electron microscope (SEM) used to characterize 
them. The entire setup (MTO and fiber-sphere) was rigidly mounted into a 
can which was hung inside a vacuum chamber by means of soft springs. A 
built-in magnètic damping vibration isolation helped in reducing the coupling 
with the main resonance of the building, occurring at 7 Hz. Finally, the vacuum 
chamber was, in turn, mounted onto an air table. This combination of vibra- 
tion isolation systems yielded peak-to-peak vibrations with Az pp < 0.02 nm 
(the detection limit in our accelerometer) for freqüències above 100 Hz. To 
reduce the damping in the MTO, the vacuum chamber was initially evacuated 
to a pressure less than 10 -5 torr. In order to reduce external vibrations, the 
chamber was then closed, the pumping station removed, and the vacuum lines 
disconnected. The pressure was kept low during the experiment by means of 
a chemical pump, made of a cold (T = 77K) activated carbon trap located 
inside the vacuum chamber. A pressure ~ 10 -5 torr was maintained with this 
setup for the length of each experimental run. When necessary, the vacuum 
chamber was re-evacuated between runs. 
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2.2. Sample characterization 



Before performing the Casimir interaction measurements, the sample was char- 
acterized using an atòmic force microscope (AFM). This serves the double 
purpose of determining the roughness of the sample and also, by repeating 
the AFM study after the Casimir measurements, to check that the sample 
was not modified during the measurement. The AFM images were obtained 
for both the plate and the sphere. The results are shown in Fig. 2. It is not 
possible for us to ensure that the region shown in Fig. 2a and Fig. 2b (before 
and after the measurements, respectively) is exactly the one giving the domi- 
nant contribution to the interaction with the sphere. We note, however, that 
we imaged with the AFM a region of (100 x 100) /iin 2 , located at the center 
of the Au-coated region of the plate, and Fig. 2a represents just a fraction of 
this region. The remainder of the investigated region, however, shows a very 
similar behavior, with identical topography except for two dust partides ob- 
served. For the sphere it is easier to find the interacting area since the fiber 
itself provides a preferential direction. Figs. 2c and 2d show AFM images of 
the bottom of the sphere, obtained before and after the experiment. As in the 
case of the plate, there are no discernible modifications in the topography. It 
is worth mentioning that a spherical surface of radius R was first subtracted 
from the data shown in Fig. 2c and Fig. 2d. For all cases shown in Fig. 2 an 
extra planarization was performed to account for the small non-linearities of 
the z-axis of the AFM's translational stage. The extra planarization for both 
the plate and sphere had \Az\ < 0.5 nm. In Fig. 3, typical cross-sections of 
the images of Figs. 2b and 2d at y = const are shown for the plate (a) and 
the sphere (b). Results obtained either before the Casimir measurements or 
at x — const are quite similar. Fig. 3 permits us to estimate the correlation 
length of surface roughness (see Section 5.4). 

These AFM images were used to characterize the roughness of the sample, 
which has to be taken into account when comparing the measured Casimir 
interaction with the theoretical predictions. As was done in [11], the roughness 
is represented by the fraction of the sample with height h^. Data resulting 
from the AFM images of Figs. 2b and 2d are presented in Figs. 4a and 4b, 
respectively. The heights hi are plotted along the vertical clXIS clS cl function 
of the fraction Wi of the total surface area having height h < hi. The width of 
each horizontal step is equal to the fraction of the total area Vi with heights 
hi < h < h i+1 (1 < i < K, where K = 105 for the plate and K — 112 
for the sphere; in both cases h\ = 0). Evidently Wi = V\ + v 2 + ■ ■ ■ + 
and w K = 1- The data of Figs. 4a and 4b will be used in Section 5.4 for 
the calculation of the roughness corrections to the Casimir pressure. As was 
noted in the Introduction, the surface roughness was considerably decreased 
compared with [11]. (The highest peaks here, hu 2 = 11.06 nm on the sphere 
and /iiQ5 = 20.65 nm on the plate, are much lower than the highest peak 
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h 59 = 98.5 nm in [11].) 



3. Experimental results for the Casimir pressure 



3.1. System calíbration 



The calibration of the system was performed in an analogous manner to that 
reported in [11]. The electrostatic force F ei (z) between the plate and the sphere 
was measured as a function of the separation z for z > 3 /im, where the Casimir 
force is smaller than 0.1% of the total force. This measurement was repeated 
for different potential differences between the plate and the sphere. For a given 
potential difference, F e i(z) between a sphere and an infinite plane is given by 



p 9^ (\r T/ ,2>fi coth(^) -ncoth(nM) 

^ sinh(nu) 



Here Va u is the voltage applied to the sphere, V is the residual potential 
difference between the metallic layers when they are both grounded, and 
cosh-u = [1 + z/R]. Fig. 5 shows the dependence of 9 (and hence the elec- 
trostatic force) on the applied voltage Va u - Adjusting Va u m Eq. (3) we found 
that F e i = ON within the experimental uncertainty for V = (17.5 ± 0.1) mV. 
Vo represents the difference in work functions between the Au layer on the 
plate and the sphere, along the path that closes the elèctric circuit between 
them. This value was observed to be constant for z in the (0.15-5) /im range, 
and it did not vary when measured over different locations in the Au layer on 
the plate. Note that F d is measured between the zero roughness levels on a 
plate and a sphere relative to which the mean vàlues of roughness profiles are 
zero (see also Section 5.4). As a result, the absolute separations from Eq. (1) 
are also determined between the zero roughness levels. For this reason, there 
is no systematic error due to roughness in the measurements of separations in 
addition to uncertainties discussed in Section 4.1 below. 

A set of 120 curves of F e i(z) was then used to fit for the parameter D = D1+D2 
(the error on D was originally estimated to be no smaller than 50 nm). To 
check the stability of the system, a subset of these curves (typically 30 of them) 
was repeated after each set of the separation dependence measurements, P(z), 
of the Casimir interaction. Three curves for different AV are shown in Fig. 6. 
The local R in the configuration when the sphere and the plane are interacting, 
and the proportionality factor k between the measured AC and 9, are also 
obtained through the measured F et (z). Finally, the following vàlues of all three 
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parameters were determined: D = (9349.7 ± 0.5) nm, R = (148.7 ± 0.2) /zm, 
and k = (50455 ± 7) N/F. The determination of D and the measurement of 9 
then yield z when z meas is measured interferometrically. 



3.2. Casimir pressure measurements 



As in [10,11], the most sensitive measurement of the Casimir interaction arises 
from measuring the change of the angular resonant frequency uj r of the oscil- 
lator in the presence of vacuum fluctuations. In this situation, u r is given by 
[10,11,18] 



b 2 OF' 
Iuj 2 dz 



(4) 



where uj q = 2ty x 702.92 Hz is the natural angular resonant frequency of the 
MTO, F is the Casimir force between a plate and a sphere, / is the moment 
of inertia of the MTO, and b 2 /I = (1.2579 ± 0.0006) ^g" 1 was obtained from 
the electrostatic interaction F e \. Using the proximity force theorem, we arrive 
at the Casimir pressure between the two plates 

P W = _^5£W. (5) 

V ; 2nR dz K ' 



The resonant frequency of the MTO at a separation z [and, hence, P{z)\ was 
measured by changing the separation between the sphere and the plate as 

z(t) — z + A z cos(cj r t), (6) 



where A z /z 1. For appropriately small vàlues of A z Eq. (4) is recovered. 
In reality, however, this is only an approximation, and the actual equation of 
motion for the MTO becomes non-linear [18] . For the determination of u r to be 
experimentally relevant, the value of A z has to be sufficiently small such that 
the effect of non-linearities is smaller than Au r , the error in the measurement 
of uj r . On the other hand, A z should be as large as possible to minimize Au r . 
The empirical approach we used in this paper is the following: At a separation 
z we determined uj r by measuring the spectral response of C r j under the 
influence of thermodynamic noise. These measurements were performed with 
an integration time r = 500 s. Subsequently a motion described by Eq. (6) was 
induced with uo = uo r fixed by means of a phase-lock-loop circuit. This was 
done while increasing the value of A z until a decrease in Au r was observed, 
while at the same time holding u r fixed at the value obtained by measuring 
the spectral response of the MTO. These results, shown in Fig. 7, allowed us 
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to define the maximum value of A z to be used at different z. The minimum 
value used was A z = 1.2 nm at z — 160 nm. 

The absolute error of u r , Au r , was found to be 2n x 6 mHz for an integration 
time of 10 s. This is approximately a factor of 1.7 smaller than in our previous 
work [11], the main reasons being a better decoupling of the apparatus from 
the environmental noise, and an improvement in the amplifier to measure C r i . 

The Casimir pressure P(z) was measured 15 times over different positions on 
the sample within a separation region z = (160 — 750) nm, and one of the 
data sets is shown in Fig. 8. Each point in the figure was obtained with an 
integration time of 10 s. 



4. Metrological analysis of experimental errors 



4-1. Random errors 



We start from the determination of the absolute error Az in the separation 
z measured between the zero roughness levels of Au films on the plate and 
on the sphere. According to Eq. (1), Az is determined by the absolute errors 
Az meas = 0.2 nm and AD = 0.5 nm. (We note from Eq. (2) that with A9 
10- 9 v^rad/Hz 1/2 and 9 < 10" 5 rad [11], the third term in the right-hand 
side of Eq. (1) makes a negligible contribution to Az.) Both z meas and D 
are assumed to be distributed uniformly within the limits of their respective 
absolute errors (note that the use of other distributions would decrease Az). 
The absolute error of the quantity obtained by the composition of the uniform 
distributions is given by [44] 



Az = min 



(Az meas + AD), k^(Az meas ) 2 + (AD) 2 



(7) 



where the correction factor kp depends on the chosen confidence level and the 
number of composed quantities (two in our case). Here and below we will use 
the confidence level (3 = 95% which results in kp — 1.1 [44]. Substituting this 
into Eq. (7) we conclude that the absolute separations are measured with an 
error Az = 0.6 nm. 

As mentioned in Section 3.2, the Casimir pressure was measured 15 times over 
the 160 nm to 750 nm separation region (with 288 to 293 points in each set 
of measurements). It is not trivial, however, to take advantage of the fifteen- 
fold repetition of the measurement in order to decrease the random error and 
to narrow the confidence interval. The difliculty is that in the measurement 
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procedure employed the separation step between any two neighboring points 
was not uniform, even in one set of measurements, and additionally was quite 
different in all sets of the data. This leads to an absence of even a few points 
(let alone fifteen) taken at the same separation which could thus be used for 
averaging in a statistical analysis. 

The usual way to deal with such data is the following: The entire separation 
range under consideration (from 160 nm to 750 nm) is divided into partial 
subintervals of length 2Az = 1.2 nm each. All points from each set of measure- 
ments are plotted together on this interval and, as a result, each subinterval 
j contains a group of a few points rij (in our case from 3 to 13). Inside each 
subinterval all points can be considered as equivalent because the value of the 
absolute separation is distributed uniformly within the limits of 2Az. 

We note that the experimental points used should be checked for the presence 
of so called "outlying" results. For this purpose inside of each subinterval j it 
is necessary to consider the quantity 

T j = —max\P j , i -P j \, (8) 

where Pj ti is the value of the Casimir pressure at point number i of the subin- 
terval j (1 < i < rij, maximum is taken with respect to i). The mean and the 
variance of the pressure are respectively defined as 

p i = ^-E s % = E( p ^ - ^) 2 - ( 9 ) 
n j i=i iL j l i=i 

If the inequality Tj > T njj i-p is satisfied, where 7^,1-/3 are tabulated quan- 
tities, then the subinterval j contains an outlying result, which should be 
rejected with a confidence probability (3 [44,45]. 

The 15 available sets of measurements were analyzed using this criterion and 
one set of measurements was found to be outlying. The subintervals containing 
points of this set satisfy the above inequality with high probability j3. As an 
example, at separations 170 nm, 174 nm, 180 nm and 250 nm the probabilities 
that the points of the rejected set are outlying are 80%, 98%, 95%, and 98%, 
respectively. For this reason in the subsequent analysis we only consider data 
from 14 sets of measurements which were confirmed not to be outlying. 

Direct calculation shows that data from these 14 measurement sets consist of 
groups (belonging to the neighboring subintervals) which are uniform in mean 
vàlues. This implies that the mean Casimir pressures Pj change smoothly in 
going from group (subinterval) j to group (subinterval) j + 1. The variance, 
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however, turned out not to be uniform in going from one subinterval to an- 
other. For this case we have utilized the theory of repeated measurements to 
implement a special procedure [46,47]. At each separation distance z , in order 
to find the uniform variance of a mean, it is necessary to consider not only 
one subinterval (which covers the value of z ) but also several neighboring 
subintervals to the left and to the right of Zo (in our case usually 4 or 5) as 
well. The specific number of such subintervals N is determined from the re- 
quirement that the variance of the mean computed over several subintervals 
be uniform in going to the next set of subintervals. 

Mathematically the confidence interval at some fixed separation z is found 
in the following way. We calculate the variance of the mean of the Casimir 
pressure in each subinterval j (1 < j < N) chosen around Zq, 



<5 p. — , 

n j\ U J V i=l '"3 



5 P 



-^U^-Pj)^^- (io) 



Then the variance of the mean at a point zq is found via the equation [46,47] 



s 2 p(zo) = max 



N 
3=1 



2„2_ 
3 S Pj 



where Xj are the influence coefficients. The maximum is taken over two sets 
of influence coefficients, Xj = 1/N, and 

Ai = -4^, (12) 

c j S Ck 
k=i 



with the constants determined from 



s 2 Pl : s 2 p 2 : . . . : Sp N = ci : c 2 : . . . : c N . 



The confidence interval at a confidence probability (3 

'Pjizo) - A rand P cxpt (z ), Pj(z ) + A rand P cxpt (^ )] (13) 

can then be found where 

A randpexpt (2o) = s _ (^) í(1+/3)/2 ( min n . _ 1). (14) 

Here t p (f) is obtained from the tabulated vàlues for the Student's t distribition 
(see, e.g., [48,49]). 
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The calculated results for A rand P cxpt (,2) over the entire interval (160-750) nm 
at /3 — 0.95 confidence are shown in Fig. 9a by the long-dashed line. As is 
seen from the figure, as the separation decreases the random error rapidly 
increases. This is a consequence of Eqs. (10) and (14) and the decrease of rij 
for subintervals j at shortest separations. 



4-2. Systematic errors 



We turn now to a discussion of systematic errors in the Casimir pressure 
measurements. By convention we call two of the errors in our experiment 
systematic which are in fact not constant, but can vary within certain limits, 
and can be described by random quantities with a uniform distribution. The 
first of these is the sphere radius with the absolute error Ai? = 0.2 jum, and the 
second is the angular resonant frequency of the oscillator u r with the absolute 
error Auj r = 2n x 6mHz (Section 3.2). 

The determination of the Casimir pressure from Eqs. (4) and (5) is in fact an 
indirect measurement. However, the error of P cxpt (z) should be expressed in 
terms of the errors of the directly measured quantities. In our case P cxpt (z) 
is the ratio of two directly measured quantities, and for this reason it is con- 
venient to use the relative errors (which we denote by 5), rather than the 
absolute errors. Eq. (7) (written in terms of the relative errors) can then be 
used once more [44] leading to 

^systpexpt^ = min L 5üür + SR ^ 1.1^(^)2+^)21 ? ( 15 ) 



where the 95% confidence level has been chosen. Recali that the quantities 
lú and b 2 /(Iul) are determined so precisely that their uncertainties do not 
contribute to the overall error. It should be noted that the first term on the 
right-hand side of Eq. (15) determines the total result at large separations of 
about 450 nm or more. At small separations S syst P cxpt is given by the second 
term on the right-hand side of Eq. (15). 

The absolute systematic error in the Casimir pressure measurements 

A systpcxpt^) = |pexpt^)| 5 systpcxpt^) ) computec l f rom Eq . ( 15 ) ; jg s h OW n in 

Fig. 9a by the short-dashed line. It is seen that at separations z > 350 nm 
the systematic error is slightly greater than the random error. At the shortest 
separations the systematic error is negligible. 
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4-3. Total experimental error 



To obtain the total experimental error in the Casimir pressure measurements, 
it is necessary to combine the random and systematic errors which are de- 
scribed by different distributions (normal or Student and uniform in our case). 
In fact, there are different methods in the literature for combining random and 
systematic errors [44]. A convenient method in practical applications [50], is 
based on consideration of the ratio 

Asyst pexpt/ \ 



At distances where the inequality r(z) < 0.8 is satisfied, the systematic error 
can be neglected and the total error of the Casimir pressure measurements at 
95% confidence level is given by 

A tot P cxpt (z) = A rand P cxpt (^). (17) 



At distances where r(z) > 8, the random error can be neglected and the total 
error at the same 95% confidence level is 

A totpexpt^ = A sy S tpexpt^_ ( lg ) 



In the intermediate region 0.8 < r(z) < 8, it is recommended [50] that one 
use the expression 



A totpexpt^ = k ^ Urandpexpt^) + ^ystpexpt^) 



(19) 



At 95% confidence level the tabulated coefficient fcg(r) varies between 0.71 
and 0.81. (For practical purposes it suffices to use hp(r) = 0.8 [50].) 

The total error on the measurements of the Casimir pressure in the separation 
region (200-750) nm is shown in Fig. 9b by the sòlid line. For the sake of 
convenience, the random and systematic errors are shown in the same figure 
by the long-dashed and short-dashed lines, respectively, on a larger scale than 
in Fig. 9a. At separations (160-200) nm the total error coincides with the 
random one, i.e. is given by the long-dashed line in Fig. 9a. As can be seen 
in Fig. 9b, at separations z > 400 nm the total experimental error is almost 
constant and varies in the interval (0.43, 0.41) mPa. 

In Fig. 10, the total relative experimental error 

fitat pexpt ^ = A tot P cxpt (2)/|P cxpt (z)| 
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is plotted versus separation. It is seen that the total relative error is almost 
constant (varies between 0.55% and 0.60%) in a wide separation region from 
170nm to 300 nm. This contrasts with previous experiments on the Casimir 
force, where the smallest relative error (1.75% at 95% confidence [22] in the 
most precise and accurate AFM experiment [6]) was achieved only at the 
shortest separation. 



5. Different theoretical approaches to the Casimir pressure be- 
tween real metals at nonzero temperature 



5.1. Lifshitz formula combined with the surface impedance or with the plasma 
model 



The dynamic determination of the Casimir pressure, as described in Section 
3.2, was carried out using a Au-coated sphere over a Au-coated plate. Since the 
thicknesses of both metal coatings were greater than the plasma wavelength 
of Au, X p = 137 nm, one can calculate the Casimir pressure as if the sphere 
and plate were composed of sòlid Au [12]. In the case of two infmite plates 
(the correction due to finite sizes of the plates is negligible, see Section 6) the 
result at temperature T is given by the Lifshitz formula [23,24,25] 



P(*) = -— É' [k±dk ±qi 



(20) 



x 



+ 



-i 



Here k± is the wave vector component in the plane of the plates, k± = \k±\, 
qf = k\ + Çf/c 2 , Çi = 2nkBTl/h are the Matsubara freqüències, ks is the 
Boltzmann constant, and riu are the reflection coeficients for two indepen- 
dent polarization states computed for the imaginary freqüències U3\ = The 
prime in Eq. (20) refers to the inclusion of a factor 1/2 for the term with / = 0. 

Eq. (20) has been derived in the framework of many different formalisms, and 
can hence be considered as the firmly established foundation for the description 
of the van der Waals and Casimir pressures between real materials with finite 
conductivity at nonzero temperature. However, the explicit expressions for 
the reflection coeflicients ry j_ are less certain. Traditionally, following Lifshitz 
[23,24,25], they were expressed in terms of the dielectric permittivity s{uj) 



h + e(i£i)qi 
h - e(iÇi)q t 



h + qi 



h - qi 



(21) 
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where kf — k\ + e{i^i)^f / c 2 . For real metals at nonzero temperature and 
separations greater than the plasma wavelength the reflection coeficients can 
also be represented in terms of the surface impedance Z [30] 

Z(i£i)cqi + 6 ] 2 (22) 
Z(^i)c qi -^\ • 1 ) 



Z(i£i)ji + cqi 
_Z(iíi)íi - cq t 



At T = the Lifshitz formula with reflection coeflicients given by Eq. (22) 
has long been known [51]. On the real frequency axis, the surface-impedance 
reflection coeflicients, Eq. (22), are commonly used in the analysis of opti- 
cal properties of metals [52]. We note that Eq. (22) contains the Leontovich 
impedance which does not depend on polarization or transverse momentum 
[31]. If the dielectric permittivity depends only on frequency (as admitted in 
the Lifshitz theory) the Leontovich impedance is given by Z{uS) = 1/ ^e{ui). 
However, the Leontovich impedance and corresponding boundary conditions 
on the surface of metal are more general since they still hold, for instance, in 
the frequency domain of the anomalous skin effect where it is impossible to 
describe the metal in terms of e{uj) due to the spatial nonuniformity of the 
electromagnètic field [52]. 

It is important to note that expressions (21) and (22) are not equivalent. They 
would only become equivalent if the Leontovich impedance Z{uS) in (22) is 
replaced by the so-called exact impedances which depend on both polarization 
and transverse momentum, 

Z\\{uM) = — U - c 2 k 2 ± , Z ± (cv,k ± ) = = , (23) 



where tv and k± are independent and not constrained by the mass-shell equa- 
tion k\ + k\ = uj 2 /c 2 vàlid for real photons [31] (k 3 is the wave vector compo- 
nent perpendicular to the plane of the plates). Within the fourth (impedance) 
approach to the thermal Casimir force it is postulated that the reflection 
properties of virtual photons on a metal boundary are the same as for real 
photons, i.e., to and k± in Eq. (23) must be constrained by the mass-shell 
equation. Under this additional condition it then follows that the reflection 
coeflicients calculated with the exact impedances (23) coincide precisely with 
(22), expressed in terms of the Leontovich impedance, at zero frequency and 
are approximately equal to (22) with a very high accuracy at all other freqüèn- 
cies (see [31] for more details). The use of reflection coeflicients (22) removes 
the contradictions with thermodynamics inherent in the first approach to the 
thermal Casimir force [26,34,35,36,37] which uses the coeflicients (21) in com- 
bination with the Drude dielectric function (see Introduction). 

Within the experimental separations (160-750) nm the characteristic angular 
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frequency of the Casimir effect £ c = c/(2z) lies within the region of infrared 
òptics. In the fourth (impedance) approach (see Introduction) the following 
calculation procedure is used to compute the vàlues of the Casimir pressure. 
First, using the tabulated optical data for Au [53], the dielectric permittivity 
along the imaginary frequency axis is computed 



*«) = l + |/í^. (24) 



Ime:(u;) is found from the tabulated complex refractive index extending for Au 
from 0.125 eV to 10000 eV [53]. To calculate e at all contributing Matsubara 
freqüències with a sumcient precision, it is necessary to extrapolate the avail- 
able tabulated data for freqüències below 0.125eV. This is usually done (see, 
e.g., [54,55]) using the imaginary part of the Drude dielectric permittivity 

eM = 1 - r^~T > ( 25 ) 



where u p is the plasma frequency and 7 is the relaxation frequency (for Au 
Up = 9.0 eV and 7 = 0.035 eV; leV = 1.51927 x 10 15 rad/s). In all compu- 
tations the following vàlues of fundamental constants, recommended by the 
National Institute of Standards and Technology, were used: c = 2.997924 x 
10 8 m/s, h = 1.05457 x 10" 34 J • s, k B = 1.38065 x 10~ 23 J/K, and e = 
1.602176 x 10- 19 C. 

Once s(i£) has been found from Eq. (24), the vàlues of the impedance for 
the imaginary Matsubara freqüències with / > 1 are calculated as Z(i£) = 

1/ \Je(iÇ). The resulting vàlues are only slightly different from those computed 
from the analytical representation for the surface impedance in the free elec- 
tron plasma model [given by Eq. (25) with 7 = 0], 

Z{%Í) = * . (26) 



The reason for minor differences is that the impedance from Eq. (26) is purely 
imaginary when evaluated on a real frequency axis. Eq. (26) ignores the small 
real part of the impedance in the region of infrared òptics arising, for instance, 
from electron-electron scattering. This real part is, of course, contained in the 
tabulated optical data. 

In the fourth approach [30,40], the value of the impedance at zero Matsubara 
frequency is obtained by extrapolation from the region of characteristic fre- 
qüències (infrared òptics in our case). As a result, the behavior Z(i£) ~ Í/u) p 
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when £ — > is obtained from Eq. (26), which from Eq. (22) leads to 



Using the resulting vàlues of Z(i£i), the Casimir pressure Pi{z) (using the 
fourth approach) was computed from Eqs. (20), (22), (24), (27) at all ex- 
perimental separations for all 14 sets of measurements at T = 300 K. As an 
illustration, several results are included in Table 1 (columns 2-4). In column 2 
the contribution of the zero-frequency term of the Lifshitz formula pf~^ (z) is 
presented. Column 3 contains the contribution of all terms Pf~ l \z) with / > 1 
in the Casimir pressure. Column 4 presents the vàlues of the total pressure 
Pa{z) (the sum of the contributions from the two previous columns). 

It is useful to compare the results, obtained by the impedance approach, with 
those from [28,29], obtained by the substitution of the plasma model dielectric 
function [Eq. (25) with 7 = 0] into the Lifshitz formula (20) with reflection 
coefficients (21) (the third approach, see Introduction). In this case at zero 
Matsubara frequency, one finds in place of Eq. (27) 



The calculated results for the contribution of the term P% (z) with l — in 
the Casimir pressure, the terms P^ l ~ 1 \z) with / > 1, and the total pressure 
Ps{z), are presented in columns 5,6,7 of Table 1, respectively. In column 8, 
for comparison, the contribution of all terms with l > 1 is calculated using the 
tabulated optical data and Lifshitz reflection coefficients, Eq. (21), expressed 
in terms of the dielectric permittivity. This column should be compared with 
column 3 where the tabulated optical data were used in combination with the 
impedance reflection coefficients from Eq. (22). 

As can be seen from Table 1 (columns 4 and 7), at the shortest separation 
z = 160 nm the plasma model gives smaller magnitudes of the Casimir pres- 
sure than the impedance approach by approximately 2.6%. The same holds 
at zero temperature [54] where the situation is quite clear and there are no 
conflicts among the different approaches. This underestimation decreases to 
only 0.7% with increasing separation until z = 250 nm. At greater separa- 
tions the differences among the data in columns 4 and 7 quickly decrease. 
The above differences at short separations are explained by the fact that the 
plasma model does not include the infernal photoelectric effect (interband 
transitions) and other processes taken into account by the tabulated optical 
data. On the other hand, if we compare columns 3 and 8 of Table 1, it is 





(27) 




(28) 
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seen that the differences are less than 0.42% (the greatest difference being at 
the shortest separation z = 160 nm). It follows that the impedance reflection 
coeficients from Eq. (22) lead to practically the same results as the Lifshitz 
coeficients given by Eq. (21). 

5.2. Lifshitz formula with zero-frequency term determined by the Drude model 
or by special prescription 

As was noted in the Introduction, the first and second alternative approaches 
to the description of the thermal Casimir force between real metals lead to pre- 
dictions which differ from those in Section 5.1. Here we compute the Casimir 
pressure in the framework of these approaches as they were formulated by 
their authors. 

We start with the first approach [26] in which the Lifshitz formula (20) and 
the reflection coeficients (21) were used to calculate the Casimir pressure 
between Au plates, taking into account the effects of nonzero temperature 
and finite conductivity. The vàlues of the dielectric permittivity at imaginary 
Matsubara freqüències e(iÇ,i) were obtained in the same way as was described 
in Section 5.1. In other words, tabulated optical data, extrapolated in the 
low frequency range by the imaginary part of the Drude dielectric function, 
were used to obtain e(i£i) by means of the dispersion relation (24). Up to 
this point the approach of [26] almost coincides with the impedance approach. 
(As was demonstrated in Section 5.1, column 8 of Table 1, the use of the 
Lifshitz reflection coeficients instead of the impedance coeficients leads to 
only minor differences at all Matsubara freqüències with l > 1.) The major 
difference in the approach in [26] is the value used for the zero-frequency term 
of the Lifshitz formula. Instead of Eq. (27), the reflection coeficients at zero 
Matsubara frequency were obtained in [26] by the substitution of the Drude 
dielectric function (25) into Eq. (21) with a result 

rf L (0,fc ± ) = l, ri >L (0,fc J _) = 0. (29) 

The approach of [26] was supported in [34,35] (see also [36,37]). 

Using the calculational procedure of [26,34,35], the Casimir pressure P\{z) 
obtained from the first approach was computed at all experimental separations 
at T = 300 K. Several results are presented in Table 2. In columns 2,3 the 
zero-frequency term P^~°\z) and the total Casimir pressure P\{z) are given, 
respectively. (Recali that the contributions of the Matsubara freqüències with 
l > 1 were already included in column 8 of Table 1.) 

We next turn our attention to the second alternative approach to the thermal 
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Casimir force suggested in [27]. This approach also used the Lifshitz formula 
(20) with the reflection coemcients (21) to calculate the Casimir pressure be- 
tween the plates made of real metals at nonzero temperature. The dielectric 
permittivity e(i£i) was found in the same way as in [26]. The major differ- 
ence between the approach [27] and those used in Section 5.1, or in [26], is a 
new alternative for the zero-frequency term of the Lifshitz formula. Instead 
of Eqs. (27), (28) or (29), it was postulated that for real metals both reflec- 
tion coemcients at zero Matsubara frequency should be the same as for ideal 
metals, i.e. 

rg(0,À; ± ) = l, rl(0,A; ± ) = l. (30) 

Using this postulate, the Casimir pressure p2(z) obtained from the second 
approach was computed from Eqs. (20), (21) at all experimental separations at 
T = 300 K. As an example, several calculated results are presented in columns 
4, 5 in Table 2, where the zero-frequency term P2~°\z) and the total Casimir 
pressure ^2(2), obtained from Eqs. (30) and (20), are given, respectively. 

Comparing column 4 of Table 1 with column 3 of Table 2, we conclude that the 
first approach leads to smaller magnitudes of the Casimir pressure than the 
impedance approach at all separations. As an illustration, the differences of the 
theoretical predictions given by these two approaches change from 2.4% of the 
Casimir pressure at z — 200 nm to 7.4% at z — 500 nm. On the contrary, if we 
compare column 4 of Table 1 with column 5 of Table 2, we conclude that the 
second approach predicts larger magnitudes for the Casimir pressure than does 
the impedance approach. The differences of the theoretical predictions given 
by these two approaches change slowly from 2.5% of the Casimir pressure 
at z — 200 nm to 2.2% at z — 500 nm. It should be emphasized that the 
calculated results in Table 2 are in agreement with those obtained by the 
authors of [26,34,35,36,37] (columns 2, 3) and [27] (columns 4, 5), where the 
first and second approaches, respectively, were developed. 

5.3. Variation of optical tabulated data and effects of spatial nonlocality 

To ensure that the theoretical approaches which use tabulated optical data 
for the complex refractive index, lead to sufnciently precise results, one should 
check that small sample-to-sample variations of the data due to, for example, 
size of the grains and impurities are not problemàtic. This was verified in 
[22] for the experiment of [6]. It was found that at z — 62 nm the possible 
variations of the optical data due to grain size lead to a 0.8% decrease of 
reflectivity, and to a decrease of less than 0.5% in the Casimir force. In our 
case the sizes of grains vary between 25 nm and 300 nm (with a mean value 
of about 150 nm), and the characteristic wavelength at z — 160 nm is À c = 
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Incj^c = 4ttz pa 2010 nm. According to [56], even for grains as small as 45 nm 
in size (the largest studied in that paper) the reflectivity at À ~ 2000 nm is 
only 0.6% less than that calculated from tabulated data. Taking into account 
that the mean size of crystallite grains in our case is several times larger than 
45 nm, we can conclude that the variations of the optical data may lead to 
a decrease in the magnitude of the Casimir pressure much lower than 0.5%, 
even at the shortest separation considered. 

One more effect to discuss, related to the finite conductivity corrections, is the 
probable contribution from the effects of spatial nonlocality (i.e., wave vec- 
tor dependence of the dielectric function). These effects influence the Casimir 
force magnitude in the region of the anomalous skin effect at large separa- 
tions, z > 2.36 /im for Au [30], a region not relevant for our experiment. An- 
other, short-separation, region where nonlocality may contribute to the van 
der Waals force, is z < \/ (Au) pa 10.9 nm which corresponds to characteristic 
freqüències £ c > u p . At such high freqüències the spatial dispersion modifies 
the freqüències of surface plasmons. In our experiment, however, the effects 
of nonlocality do not contribute to the Casimir pressure. This can be un- 
derstood by noting that the largest characteristic frequency here, calculated 
at z — 160 nm, is £ c = 9.4 x 10 14 rad/s, which is almost 15 times smaller 
than uj p = 1.37 x 10 16 rad/s. A direct calculation shows that the Matsub- 
ara freqüències Ço, £i, . . . ,£35 contribute 99.96% of the total Casimir pressure 
at z — 160 nm, and even £ 35 = 8.61 x 10 15 rad/s < u p . We note in passing 
that the contribution of the nonlocal effects for Au of about 2% at a sepa- 
ration z = \ p = 137 nm, claimed in [57], is in error, as was recognized by 
the authors in a recent Erratum [58]. (In fact the formalism of this reference 
leads to a correction of only 0.2% due to spatial nonlocality at z — X p ; the 
confusion arises from a notation L/\ p on the rr-axis in Fig. 2 of [57] instead 
of 2nL/\ p .) Furthermore, in [57] a spatially nonlocal dielectric function was 
used in the frequency region of infrared òptics where it is in fact local [52,59]. 
As was specifically demonstrated in detail in [60], the short-separation effects 
of spatial nonlocality connected with the surface-plasmon charge fluctuations 
contribute considerably to the force at separations of the order of vf/oj p ~ 1 A, 
where vf is the Fermi velocity. It was explicitly proved in [60] that for much 
larger distances z ^> Vp/üJp (in fact for z larger than several nanometers) 
the usual local Lifshitz formula is applicable. This is in agreement with our 
conclusion that at separations z > 160 nm there is no noticeable contribution 
from the effects of nonlocality to the Casimir pressure. 



5.4- Casimir pressures taking account of surface roughness 

At separations below 1 /im surface roughness corrections may contribute from 
a fraction of a percent up to 20% of the measured force in different experi- 
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ments [12]. For this reason they should be carefully taken into account in pre- 
cision computations of the Casimir pressure. The contribution of roughness 
is calculated with the help of the AFM characterization of both surfaces (see 
Section 2.2). As is seen in Figs. 2 and 3, the roughness is mostly represented 
by stochastically distributed distortions with a characteristic lateral size of ap- 
proximately / ~ (500 — 600) nm. This size is larger by a factor of ps 3 than z at 
the shortest separations (160-200) nm, where the roughness corrections may 
give a noticeable contribution to the Casimir pressure in our experiment. The 
influence of roughness can be calculated by the additive method of geometrical 
averaging [12]. As shown in [61], under the condition z < l the additive method 
gives approximately the same results as more detailed methods. Moreover, the 
correction due to nonadditivity can be simply estimated. 

There are two modifications to the additive method. The first one is based 
on the use of the Casimir pressures Pi(z), found in the above four theoretical 
approaches as calculated in Sections 5.1, 5.2. Let us first determine the zero 
roughness levels H^ ,s on a plate and a sphere, respectively, relative to which 
the mean vàlues of the functions, describing roughness, are zero (note that 
the experimental separations are measured between the zero roughness levels 
[12]): 

f; (#r - = o. (si) 
i=i 



Here the quantities vf ,s and hf s are taken from Figs. 4a and 4b for the plate 
and the sphere, respectively (K was defined in Section 2.2). To find the Casimir 
pressure taking into account all three corrections, the vàlues P(z) should be 
geometrically averaged over all possible separations between the rough sur- 
faces, weighted with the probability of each separation [3,11,12] 

105 112 

p ^ or{z) = j2^ viv'Pt (z + HI + H s - h\ - hi) . (32) 

fc=l j=l 



Note that this expression does not reduce to a simple multiplication of the 
correction factors due to nonzero temperature and finite conductivity on the 
one hand, and due to surface roughness on the other, but takes into account 
their combined (nonmultiplicative) effect. 

The second one is a simpler additive method. It was used for calculating 
the corrections due to stochastic surface roughness in [62]. According to this 
method, the Casimir pressure between two ideal metal parallel plates covered 
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by stochastic roughness with variances 8 p st and 8 s st is given by 
P r (z) = rj r P (z), r] r = l + 10 



z 



z 



(33) 



where P (z) = —TT 2 hc/ (240^ 4 ) is the Casimir pressure between perfectly shaped 
parallel plates. The variances describing the stochastic roughness are found 
from the formula 



K 



mr = e ( HP o' s - hrr 



(34) 



í=i 



Using data from Figs. 4a, 4b, one obtains the vàlues for variances 5 p st = 4.06 nm 
and 8g t = 1.91 nm. The smallness of these variances justifies the neglect of 
the fourth order terms on the right-hand side of Eq. (33) (the fourth order 
contributes less than 0.01% even at the shortest separation z = 160 nm). We 
can now obtain the complete theoretical result, including the corrections due 
to surface roughness, nonzero temperature and finite conductivity, by the use 
of a multiplicative procedure 

P^ \z)=r lr P l {z). (35) 



It is interesting to compare the theoretical results for the Casimir pressure 
taking into account all corrections contained in Eq. (32), and the simpler ex- 
pression in Eq. (35). As an example, at the shortest separation z = 160 nm, 
Eq. (32) gives P 4 thcor = -1.1515 Pa and Eq. (35) results in P 4 thcor = -1.1531 Pa, 
a difference only 0.14% of the Casimir pressure. At z — 200 nm the complete 
Casimir pressures are P^ heor = —0.51143 Pa and p* heor = —0.51185 Pa com- 
puted from Eqs. (32) and (35), respectively. In this case the difference is even 
smaller: 0.08% of the Casimir pressure. For larger separations the role of the 
nonmultiplicative effects decreases and one may use Eq. (35) to calculate the 
theoretical vàlues of the Casimir pressure with sufficiently small uncertainty. 
The same relative sizes of roughness corrections are obtained for the other 
theoretical approaches. 

Note that the role of the surface roughness in the improved experiment is 
rather modest. Using the data of Tables 1 and 2 we find that at z — 160 nm 
roughness contributes 0.65%, and at a separation z = 200 nm only 0.42% of 
the Casimir pressure. 

It is important to bear in mind that both Eqs. (32) and (35) are based on the 
approximation and do not take into account diffraction-type effects arising in 
the case of periòdic roughness with small periods l < z [61], or stochastic 
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roughness with small correlation length [63]. To obtain an estimate for the 
upper limit of the contribution of diffraction-type effects we follow [22]. For 
this purpose, the bound on the roughness correlation length is found to be 
horr > 600 nm, where 600 nm is the lateral size of the largest cluster on the 
surface. We next consider periòdic roughness with a period l corr = 600 nm 
using the fact that diffraction-type effects are greater for a periòdic function 
with a period l corr than for a random function with a correlation length l corr . 
The diffraction-type effects for a periòdic function can be computed in the 
framework of the functional approach of [61]. At separations z = 160 nm and 
300 nm one obtains z/l corr = 0.27 and 0.50, respectively. The coefficient c corr 
in the expression 



tf orr = 1 + 10c c 




(36) 



taking diffraction-type effects into account, is related to that plotted in Fig. 2 
of [61] by the equation 

r - r - 1 ~ dCcorr ("VT) 

^corr ^corr -° ■ V / 



(In [61] the results for the Casimir energy between plates are presented, not 
those for the pressure.) Using the data of Fig. 2 from paper [61] and Eq. (37), 
we find c corr ~ 1 and c corr fa 1.18 at separations z = 160 nm and z = 300 nm, 
respectively. Hence at z — 160 nm there is no noticeable contribution from 
correlation effects to the roughness correction. From Eqs. (33) and (36) it 
follows that at z = 300 nm rj r = 1.0022 and r]™"' = 1.0026, leading to a 
contribution from correlation effects to the roughness correction of less than 
0.04%. Because of this, the correlation effects of roughness in the experiment 
under consideration can be ignored. 



5. 5. Contribution of patch potentials 



In the experimental configuration described in Sections 2, 3, the electrostatic 
force due to the residual potential difference between the plate and the sphere 
is negligible. It is conceivable, however, that spatial variations of the surface 
potentials due to grains of polycrystalline metal film (the so called "patch 
potentials") simulate the Casimir force [64]. Here we use the general results 
of [64] in order to demonstrate that the patch effect is negligible in our exper- 
iment. According to [64], for a configuration of two parallel plates the elèctric 
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pressure due to random variations in patch potentials is given by 



2e Q a 2 v 



k, 





2 - kl : 



/ 



dk. 



(38) 



k. 



Here <t„ is the variance of the potential distribution, and k max (& mm ) are the 
magnitudes of the extremal wavevectors corresponding to minimal (maximal) 
sizes of grains. Using the work functions of Au for different crystallographic 
surface orientations (Vi = 5.47 eV, Vi = 5.37eV, and V3 = 5.31 eV), and 
assuming equal areas of respective crystallographic planes, we obtain 



As was mentioned in Section 5.3, the extremal sizes of grains in the Au layers 
covering the test bodies are \ m in ~ 25 nm and X max — 300 nm. This results 
in k max = 0.251 nm -1 and k m i n = 0.0209 nm -1 . Substituting these vàlues into 
Eq. (38) yields the additional pressures due to the patch elèctric forces. At 
the shortest separations z = 160 nm and z = 170 nm the patch pressures are 
ppatch _ 0.42 mPa and 0.25 mPa, respectively. In comparison to the Casimir 
pressures at the same separations, the relative contributions of the patch effect 
to the pressure are 0.037% and 0.027%, respectively, and further decrease with 
increasing z. It follows that patch effects do not play any significant role in 
our Casimir pressure measurements by means of a micromachined oscillator. 



6. Metrological analysis of theoretical errors 

As it is evident from the foregoing, the theoretical computation of the Casimir 
pressure with all relevant corrections is a rather complicated procedure requir- 
ing a variety of data. Bearing in mind the eventual comparison of theory with 
experiment, it is important to analyze all uncertainties introduced at different 
stages of the computations. 

According to Eq. (5), the derivative of the Casimir force between a plate and 
a sphere with respect to separation is related to the Casimir pressure between 
two plane plates by the use of the proximity force theorem. It is well known 
that this "theorem" is approximate and its relative error S p f t P theor is less than 
z/R [65,66,67]. This error is separation-dependent and increases from 0.107% 
at a separation z = 160 nm to 0.5% at z — 750 nm (see the short-dashed line 
in Fig. 11). 




i=l 



(39) 
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The computation of the Casimir pressure taking into account nonzero tem- 
perature and finite conductivity in Sections 5.1, 5.2 was based on the use of 
tabulated optical data for Au (approaches 1, 2 and 4). The primary source 
of the theoretical errors is the sample to sample variation of these data. As 
was shown in Section 5.3, in the experiment under consideration the variation 
of the optical data leads to an uncertainty in the magnitude of the Casimir 
pressure which is much smaller than 0.5%. To be conservative, we admit an 
uncertainty in the computation of the finite conductivity corrections as large 
as 5 c P· thcor = 0.5% over the entire measurement range. In fact, this error also 
includes all differences arising when one uses the Lifshitz reflection coeflicients 
(21) instead of impedance coefficients (22) and vice versa (see Section 5.1). 
The theoretical error 5 C is shown in Fig. 11 by the long-dashed line. 

The other uncertainties discussed above, such as the contribution of nonlo- 
cality to the effect of finite conductivity (Section 5.3), the diffraction-type 
contributions to the effect of surface roughness (Section 5.4), and the correc- 
tion to the Casimir pressure due to the patch potentials (Section 5.5), were 
shown to be far smaller than the previous two discussed above, <5 p /tP ! theor and 
( s > c P/ heor . For this reason they can all be neglected. One further correction to 
the Casimir pressure, arising from the finite size of plates, is also negligible. 
As was shown in [11], in the plate-sphere configuration of a micromachined 
oscillator this correction is less than 0.04% at z = 500 nm, becoming smaller 
as the separation decreases. But even this correction does not contribute in 
our case of two equivalent parallel plates because its derivative with respect to 
separation vanishes. In fact, to obtain the correction due to the finiteness of 
the plates, one should consider corrections of higher order in z/R for a sphere 
above a plate. As a result, at least an extra power of the small parameter 
z/R appears in the correction, making it completely insignificant in the error 
analysis. 

Let us now determine the theoretical error resulting from combining the two 
major errors <5 p / t Pj theor and <5 c P í thcor . Both random quantities are described by 
the same distribution law which is close to a uniform distribution. For this 
reason the method of [44], already used in Eqs. (7), (15), can be applied once 
more, giving 



5 P^(z) = min [(8 pft PÏ bBOt + <W heor ) , 



(40) 



As in Eqs. (7) and (15), the 95% confidence level is chosen. Note that in this 
case the minimum is achieved using the second term on the right-hand side of 
Eq. (40). 



Bearing in mind the ensuing comparison of theory with experiment, it is perti- 
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nent to consider the error in the theoretical pressures resulting from the exper- 
imental error in the determination of separation distances Az [68]. Although 
this error is not purely theoretical (theoretical pressures by themselves can be 
calculated at any exact separation distance), it is of great importance when 
one computes the theoretical pressure at an experimental point defined with 
an error Az. Using the main theoretical dependence of the Casimir pressure 
on the inverse fourth power of the separation one obtains 5 z Pj theor = AAz/z, 
where Az = 0.6 nm at 95% confidence [see Eq. (7)]. The value of 5 z P* heor 
changes from 1.5% at z — 160 nm to 0.32% at z — 750 nm (see the dotted line 
in Fig. 11). Thus, this error is the primary one (especially at short separations) 
arising in the comparison of theory with experiment. 



In [22] an additional fitting procedure was proposed in order to decrease the 
error in the determination of absolute separations down to Az = 0.15 nm. 
For this purpose, for each set of measurements, z-values of all points were 
simultaneously changed within the limits of their joint absolute error in order 
to minimize the root-mean-square deviation between theory and experiment. 
This procedure, which is quite reasonable at very short separations (z > 60 nm 
in [22]), is not applied in our conservative metrological analysis because the 
method of least squares is not well-adapted for cases when both the argument 
(separation distance) and function (pressure) are determined with comparable 
relative errors. 



We are now in a position to determine the total theoretical error of the Casimir 
pressure computations. To do this, we combine the errors determined above, 
5 Pï hcor and <5 2 Pf cor . To use the analogy with Section 4.3, the error S P^ heoT 
can be considered as "random" (in the sense that it is not described by a 
uniform distribution), and the error S z P^ heoT can be likened to "systematic" 
(as described by a uniform distribution). We then find ourselves in the regime 
of applicability of Eq. (19), and at 95% confidence, we obtain the final result 



çtot^theor^ = Q g 5 z P^° V {z) + 5 P^ C0I (z) . (41) 



This result is represented by the sòlid line in Fig. 11. It is seen that the total 
theoretical error decreases from 1.65% to 1% when the separation increases 
from 160 nm to 380 nm. With further increase of separation the total theoret- 
ical error decreases more slowly from 1% to 0.9%. 
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7. Comparison of theory and experiment and stronger constraints 
on thermal effects 



7.1. Measure of agreement between theory and experiment 



In all previous experiments on the Casimir effect [2,3,4,5,6,7,8,9,10,11,12] the 
agreement between theory and experiment was discussed in terms of the root- 
mean-square deviation between experimental and theoretical vàlues of the 
force over the entire range of separation distances. In [7] this approach was 
criticized as inadequate when the force decreases rapidly as the separation 
distance increases, although no alternative approach was suggested. Here we 
compare experiment and theory with the help of rigorous metrological criteria, 
not based on the concept of root-mean-square deviation. 

From Sections 4.3 and 6, we have at our disposal two independently obtained 
confidence intervals at the same confidence probability of 95% (the first for 
experiment and the second for theory). Both confidence intervals were deter- 
mined by the application of a common statistical procedure, and the resulting 
errors for both experiment and theory can be described by random variables 
with a common distribution law. With this in mind, we consider the new 
variable P£ hcoT (z) — P cxpt (z) and determine the total absolute error on this 
quantity, and hence the confidence interval at a confidence probability 95%, 
using the same composition law as in Eqs. (7) and (15), i.e. 



\tot 



PÏ hcor (z) - P cxpt (z)] = min { [A tot P cxpt (z) + A tot i^ thcor (z)] , (42) 



l.l\[A tot P expt (z)Y + A tot Pl hem (z) 



Here A tot P cxpt (,2) is given by the sòlid line in Fig. 9b, and A tot Pj theor (z) = 
pthcor ^totp-thcor^ where ^totp-thcor ig g^own by the sòlid line in Fig. 11. Note that 

the total absolute theoretical error is almost independent of the theoretical 
approach. For this reason we omit the index i = 1, 2, 3, 4 in the notation 
for a confidence interval. The minimum here is achieved using the second 
contribution from the right-hand side of Eq. (42). We also point out that to 
be conservative we use the value kp — 1.1, as would describe the composition 
of two uniform distributions (for other distributions the value of kp can be 
somewhat smaller). 

The confidence interval for the quantity [Pl hcor (z) — P cxpt (z)} at 95% confidence 
probability is given by 



^tot ^pthcor ^ _ pvxpt^z^j ^tot ^pthcor _ pcxpt^ 



(43) 
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The meaning of the confidence interval (43) is that, if the theory is in agree- 
ment with experimental data, the mean value of [Pl heor (z) — P cxpt (z)] must 
belong to this interval with a 95% probability. This criterion can be used for 
comparison of the above four theoretical approaches with experiment. 



7.2. Comparison of dífferent theoretical approaches with experiment 



We start from the approach based on the surface impedance (the fourth ap- 
proach). The sòlid lines in Fig. 12 exhibit the confidence interval (43) obtained 
from Eq. (42) within the separation ranges (160-250)nm (a) and (250-750)nm 
(b). In the same figure the differences of the theoretical (based on the surface 
impedance) and experimental Casimir pressures are plotted for one typical 
set of 14 measurements. It is evident that the true value of the quantity un- 
der consideration p^ heor — p ex pt (which can only be achieved for a complete 
correct theory and an infinite number of measurements) is zero. As an exam- 
ple, for the measurement set plotted in Fig. 12, (pt hcor — pcxpt^ _ 0.008 mPa, 
where the averaging was taken over all 288 points belonging to this set. If 
only points at separations z > 250 nm are considered (242 points), one ob- 
tains (p* hcor — p cx pt^ — 0.09 mPa. These vàlues should be compared with 
the half-width of the confidence interval which decreases from 27.7 mPa at 
z = 160 nm to 0.46 mPa at z — 750 nm. Remarkably, only 9 of 288 points (i.e. 
3.1%) fall outside the confidence interval in Fig. 12. 

In Fig. 13a,b the same confidence interval (43) is shown once more and the 
points p* heor — p cx pt are plotted utilizing all 14 sets of measurements. For 
all 4066 points in Fig. 13 one has (P4 heor — P^p*) = 0.11 mPa. If only separa- 
tions z > 250 nm are considered (3442 points) the result is (P4 heor — p ex P^ = 
0.12 mPa. Although these vàlues are less than the half-width of the confidence 
interval, they are not sufficiently informative. In the case of Fig. 13, where all 
points are plotted, it is much more important to consider the local vàlues of the 
quantity APj(z) = {Pl hcor (z) — P cxpt (z)), where the averaging is done over the 
vicinity of a point z, in comparison with the local half-width of the confidence 
interval. Performing the averaging over the intervals (zí — 3 nm, Zi + 3 nm) (so 
that each contains approximately 43 points) at different z = Zi one obtains the 
vàlues for the mean differences between the theoretical (the fourth approach) 
and experimental Casimir pressures included in column 3 of Table 3. For com- 
parison, the vàlues of the half-width of a confidence interval are included in 
column 2. Columns 4-6 are discussed below. 

Comparing columns 2 and 3 of Table 3 we notice that all mean vàlues of 
the differences between theoretical and experimental Casimir pressures are 
well within the half-width of the confidence interval. The number of separate 
points which fall outside of the confidence interval in Fig. 13 is 207, i.e. 5.09% 
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of the total number of points from the fourteen sets of measurements. This is 
in accordance with expectations from the value of the confidence probability. 
To conclude, the experimental data of Section 3.2 are proved to be in excel·lent 
agreement with the theoretical computations of the Casimir pressure in the 
framework of the impedance approach (see also [41,42] arriving to the result 
in agreement with this approach). 

Quantitatively the agreement between theory and experiment can be charac- 
terized by the quantity 



A 1 



tot 



* t0t W = 1 lnthror/ „ L - (44) 



|ptheor^)| 



Importantly, in the experiment under consideration 5 tot (z) depends only slightly 
on z in a wide separation range, decreasing (with the fourth theoretical ap- 
proach) from 5 tot {z) = 1.9% at z = 170 nm to 5 tot {z) = 1.4% within the 
270 nm < z < 370 nm interval, and then increasing to 5 tot (z) = 1.8% at 
z = 420 nm. The largest vàlues of <5 tot at z < 170nm and z > 420 nm are 
<5 tot (^) = 2.4% at z = 160 nm and 5 tot (z) = 13% at z = 750 nm. Thus, the 
Casimir pressure measurement by means of a micromachined oscillator is the 
first Casimir effect experiment where metrological agreement between theory 
and experiment at the 1.5% level has been achieved at 95% confidence in a 
wide separation region. 

The results of a comparison between experiment and the third theoretical 
approach, based on the plasma dielectric function, are presented in column 4 
of Table 3. Comparing column 4 with column 2, one can conclude that the 
third approach is also consistent with the experimental data because all mean 
vàlues of the Casimir pressures belong to the confidence interval. At the same 
time, it is seen that at the shortest separations (170, 250) nm this approach 
does not agree with experiment as well as does the impedance approach. The 
reasons for this were discussed in Section 5.1 (see also [54]). 

We next compare the first theoretical approach with the experiment. In Fig. 14a, 
we plot the differences of the Casimir pressures P* heor — p ex pt versus separa- 
tion z (where P* heor was computed in the approach of [26,34,35] as explained 
above) for all 14 sets of measurements (i.e., in the same way as in Fig. 13 
for the surface impedance approach). The sòlid line gives the half-width of 
the confidence interval A tot (pt heOT - P cx p*) « A tot (P 4 thcor - P cxpt ) at 95% 
confidence (note that there are practically no points below the z-axis). As 
is seen from Fig. 14a, in a wide separation region 230 nm < z < 500 nm all 
points p^ heo1 — p ex pt f a N outside the confidence interval. It follows that the 
theoretical prediction of the first approach is excluded by experiment at 95% 
confidence. 
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In column 5 of Table 3 the mean vàlues of the differences between theoretical 
(obtained from the first approach) and experimental Casimir pressures are 
presented at different separations from 170 nm to 700 nm using the data from 
all 14 sets of measurements. The comparison with the half-width of the confi- 
dence interval (column 2) shows that at all separations the mean differences of 
the Casimir pressures fall outside the confidence interval. It follows that the 
first theoretical approach is excluded by the experiment at 95% confidence 
level within an even wider separation region, 170 nm < z < 700 nm. 

In fact, the above confidence interval (43) was obtained in a rather conserva- 
tive manner. The comparison of data from columns 2 and 5 of Table 3 shows 
that even if the confidence interval were to be widened to achieve 99% confi- 
dence probability, the quantities (p* heor — pexpt^ wou i c [ s ^[\\ rema í n outside this 
interval within some separation region. To make this argument quantitative 
we calculate the half-width of a new confidence interval from the equality 

Atot l pthcor pcxpt\ , / \ 

^0.99 ^1 f ) _ í[l+0.99]/2(2) _ 2 31 ^ 



A 0°95 (P! the ° r - P cxpt ) í[l+0.95]/2(2) 



(Being conservative, we preserve only two degrees of freedom, i.e. the minimum 
value from the analysis of Section 4.1). Using the data from column 2 in Table 3 
for = A tot , one obtains (pf ™ - P cxpt ) = 3.67 mPa, 1.46 mPa, 

and 1.13mPa at separations z = 300 nm, 400 nm, and 500 nm, respectively. 
Comparing these results with the vàlues of (^p^ hcor — p cx P^ from column 5 in the 
same table, one concludes that the quantities (P 1 theor — pcxpt^ g ^jj £ a jj ou t s ide 
the new confidence interval in the separation region 300 nm < z < 500 nm. 
We conclude that the first theoretical approach to the thermal Casimir force 
is excluded experimentally with 99% confidence. 

Finally, we compare the second theoretical approach [27] with the experi- 
mental data. In Fig. 14b, we plot the differences of the Casimir pressures 
ptheor _ pexpt versus separation z for all 14 sets of measurements, as in 
igs. 13 and 14a. The sòlid lines show the half-width of the confidence in- 
terval A tot (P 2 thcor - P ex Pt) ~ A tot (p 4 theor - p cx pt) at 95% confidence (there 
are no points above the z-axis in the interval 160 nm < z < 420 nm in 
the approach of [27]). As is seen from Fig. 14b, within the separation re- 



pthcor pexpt 



fall outside the 



gion 160 nm < z < 350 nm almost all points 
confidence interval. It follows that the second theoretical approach is excluded 
experimentally at 95% confidence. 

This conclusion is confirmed by considering the mean vàlues (P2 hcor — P^p*) 
obtained from the results from all 14 sets of measurements. In column 6 of 
Table 3 the mean differences (pt heor _p ex P^ are presented at different separa- 
tions from 170 nm to 700 nm using the data from all 14 sets of measurements. 
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The comparison of these data with the half-width of the confidence interval 
in column 2 shows that at separations 170 nm < z < 350 nm the mean dif- 
ferences of the theoretical and experimental Casimir pressures fall outside the 
confidence interval. This confirms the conclusion that the second theoretical 
approach is excluded by the results of the Casimir pressure measurements 
using a micromachined oscillator at 95% confidence. 



8. Stronger constraints on long-range interactions 

As mentioned in the Introduction, measurements of the Casimir force between 
metals have been successfully used to obtain stronger constraints on hypo- 
thetical long-range interactions. Such interactions have long been predicted 
in elementary particle physics (see literature cited in [11]). They can arise 
from the exchange of light elementary partides (scalar axions, graviphotons, 
dilatons, and moduli among others [69,70]) predicted by many extensions of 
the Standard Model, and as a consequence of extra-dimensional theories with 
low compactification scales [71]. In most cases the potential energy between 
two point masses m\ and m 2 separated by a distance r is given by the usual 
Newtonian potential with a Yukawa correction 



where G is the gravitational constant, a is a dimensionless constant charac- 
terizing the strength of the Yukawa force, and À is its interaction range. 

There have been many proposals to constrain a and À over sub-/xm inter- 
action ranges from van der Waals and Casimir force measurements between 
macroscopic bodies (see references on the subject in [72]). During the last 
few years the previously known constraints on (et, À) were strengthened by 
up to 4500 times from modern Casimir force measurements between metals 
[13,14,15,16,17]. In [11] an additional strengthening by up to a factor of 11 
was achieved from the first precise Casimir pressure measurement by means 
of a micromachined oscillator [10,11]. 

The results presented above from improved measurements of the Casimir pres- 
sure compared with different theoretical approaches allow us to improve the 
constraints found in [11]. Previous limits have been strengthened by up to a 
factor of 2, and their application range widened, while increasing their reliabil- 
ity This has been made possible by the use of the above metrological analysis 
which confirmed the agreement between experiment and two traditional the- 
oretical approaches within a wide separation range. (In previous experiments 
constraints were usually obtained at some fixed separation distance.) 




(46) 
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As in [11], the Newtonian gravitational force between the plate and the sphere 
of a micromachined oscillator, and the equivalent gravitational pressure be- 
tween two plates are negligible. To calculate the equivalent Yukawa pres- 
sure between two plates, one needs the detailed structure of the sphere and 
plate materials. In the experiment under consideration an AI2O3 sphere of 
density paio — 4.1 x 10 3 kg/m 3 was coated with a layer of Ti of thickness 
Atí = 10 nm with ptí = 4.51 x 10 3 kg/m 3 , and a layer of Au of thick- 
ness A s Au = 200 nm with p Au = 19.28 x 10 3 kg/m 3 . The Si plate of den- 
sity Psí = 2.33 x 10 3 kg/m 3 was first coated with a layer of Pt of thickness 
A Pt = 10 nm, with p Pt = 21.47 x 10 3 kg/m 3 , and then with a layer of Au of 
thickness A p Au = 150 nm. Considering that the conditions z, X <C R, L are 
satisfied (where L = 3.5 /im is the thickness of the plate), the equivalent hy- 
pothetical pressure between two parallel plates with the above layer structure 
is [11] 

P h yp( z ) = -2TxGa\ 2 e' z ' x 

x pau ~ (PAu - pn) e-^' x - (p Ti - Paio) e-( A ^+ A «)A] (47) 
x pau - (Pau - PPt) e~ A - /A - (p Pt - psi) e'^ +A ^ x ] . 

Note that surface roughness, which was substantially reduced in this experi- 
ment, cannot significantly affect the magnitude of a hypothetical force with 
an interaction range of about 100 nm. 

We can now obtain constraints on the hypothetical Yukawa pressure from the 
agreement between our measurements of the Casimir pressure and theory at 
95% confidence. According to our results, no deviations between traditional 
theories and experiment were observed, i.e., the hypothetical pressure should 
be less than or equal to the half-width of the confidence interval 

|P hyp (^)| < A tot \P thcor {z) - P cxpt {z)} . (48) 



Note that Eq. (48) is metrologically reliable, since it does not use the root- 
mean-square deviation as a measure of agreement between theory and exper- 
iment to obtain constraints on hypothetical interactions. It is also worthwhile 
to note that the Yukawa-type hypothetical interaction depends on the sep- 
aration distance quite differently than the thermal Casimir pressures of the 
first and second approaches, which were excluded experimentally in Section 
7.2. Because of this the mimicry of one phenomenon by another is extremely 
unlikely. 

The numerical analysis of Eqs. (47) and (48) leads to the conclusion that the 
strongest constraints are obtained within the interaction region 40 nm < A < 
370 nm from the measurement data at separations from z = 210 nm (where 
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the half-width of the confidence interval is 6.89 mPa) to z — 450 nm (where the 
half-width is equal to 0.53 mPa). In fact, with the increase of À, the strongest 
constraints on a are obtained from the measurement data at larger separa- 
tions. The resulting constraints on a are plotted in Fig. 15 for different vàlues 
of À (line la). In the same figure constraints from earlier experiments are also 
shown. These were obtained from previous experiments with a micromachined 
oscillator (line lb), from old measurements of the Casimir force between di- 
electrics [12] (line 2), from Casimir force measurements by means of a torsion 
pendulum [2,13] (line 3), and by the use of an atòmic force microscope [6,16] 
(line 4). In all cases the region in the (a, A) plane above the line is excluded, 
and below the line is allowed by the experimental results. Note that the con- 
straints from our experiment, line la, are found at 95% confidence. For the 
previous constraints, given by lines lb, 2, 3, and 4, the confidence levels were 
not determined. 

As is seen from Fig. 15, the present experiment leads to the strongest con- 
straints in the interaction range 40 nm < A < 370 nm which is wider than 
in the experiment of [11]. Comparing the constraints from lines 2 and 4, the 
largest improvement (by a factor of 20) is achieved at A ps 150 nm. If we com- 
pare the present results with those of [11], the largest improvement, by a factor 
of 2.2, is achieved at A ~ 316 nm. It should be noted that for the first time 
the new constraints, given by line la, completely fill in the gap between the 
modern constraints obtained by the atòmic force microscope (line 4) and those 
obtained using a torsion pendulum (line 3). Furthermore, as was noted above, 
the constraints of line la are determined at the 95% confidence level, which 
makes them the most reliable constraints obtained to date from measurements 
of the Casimir force. 



9. Conclusions and discussion 

In this paper we have presented the results of an improved experiment on the 
dynamical determination of the Casimir pressure between two plane parallel 
Au coated plates using a micromachined oscillator. (We recali that the pres- 
sure between two plane plates was inferred from the Casimir force between a 
plate and a sphere using the proximity force theorem.) Many improvements in 
the previously performed experiment of [11] were made yielding a large divi- 
dend in precision, accuracy and reliability of the results. The chief advantages 
of the new experiment lie in a great improvement of the surface roughness, and 
a decrease by a factor of 1.7 in the error in determining of the absolute sepa- 
rations. The metrological analysis of all experimental errors was performed at 
95% confidence. In doing so, the data were analyzed for the presence of out- 
lying measurements. This permitted us to find the total experimental error as 
a function of separation for subsequent comparison of experiment and theory. 
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To utilize all the advantages of the improved measurements, the Casimir pres- 
sures between two Au plates were calculated in the framework of four theoreti- 
cal approaches proposed in literature based on the Lifshitz formula at nonzero 
temperature, with the reflection coefficients expressed in terms of both the 
surface impedance and the dielectric function. The complete optical data for 
Au were utilized to take into account both the finite conductivity and the 
thermal corrections to the Casimir force. In doing so, many relevant factors 
and properties of Au films were computed or estimated, such as the variation 
of the tabulated optical data due to the grain structure of a metal, and the 
influence of spatial nonlocality and patch potentials. The surface roughness 
was carefully taken into account including the role of nonmultiplicative and 
diffraction-type effects. The metrological analysis of all theoretical errors was 
performed at 95% confidence (including the error resulting from use of the 
proximity force theorem, finiteness of the plate area, and errors in the deter- 
mination of the surface separations). As a result, the total theoretical error 
(independent of the experimental one) was combined using the metrological 
criteria. 



To compare different theoretical approaches and experiment, we replaced the 
root-mean-square deviation previously used in literature, by a rigorous pro- 
cedure for the composition of the experimental and theoretical errors at 95% 
confidence. As a result, a confidence interval was found for the random vari- 
ables Pf cieor {z) — P oxpt (z) , where i = 1, 2, 3 and 4 denotes the different the- 
oretical approaches described in Sections 1,5. It was demonstrated that with 
the impedance approach and the plasma model the mean vàlues (P£ h £° r (z) — 
P cxpt (z)), computed at all separations, fall inside the confidence interval. As 
a consequence, these approaches are consistent with experiment at 95% confi- 
dence. In terms of the relative uncertainty, there is agreement between each of 
them and experiment at the level (1.4-1.9)% when the separation changes in 
the interval 170 nm < z < 420 nm. This is the best agreement achieved to date 
at 95% confidence compared to any other experiment measuring the Casimir 
force. Despite the fact that the thermal Casimir force computed within the 
impedance approach or the plasma model was found to be consistent with ex- 
periment, the thermal correction by itself has not yet been directly measured. 
According to these approaches, at T = 300 K this correction is small, in qual- 
itative agreement with case of ideal metals, and can be readily measured in 
the near future by means of proposed experiments [73,74]. 

The experimental results were used as a test for the alternative approaches to 
the thermal Casimir force which predict large thermal corrections at short sep- 
arations. The above metrological comparison between experiment and theory 
was repeated with the other vàlues Pl hcor (z), computed in accordance with the 
first approach [26,34,35]. It was found that at all separations from z = 170 nm 
to z = 700 nm the mean vàlues (Pl heoT (z) — P cxpt (z)) fall outside the con- 
fidence interval. Moreover, at separations 300 nm < z < 500 nm the mean 



36 



Casimir pressures computed in the first approach fall outside the even wider 
99% confidence interval. This enabled us to conclude that the first theoretical 
approach to the thermal Casimir force is excluded by our measurements at 
99% confidence. 

The same metrological comparison was carried out with respect to the second 
approach [27]. The Casimir pressures P^^i^z) were computed and compared 
with experiment. It was shown that the mean vàlues {P2 hcor (z) — P cxpt (z)) 
fall outside the confidence interval for separations 170 nm < z < 350 nm. 
Thus, the second theoretical approach is also excluded experimentally at 95% 
confidence. 

Finally, the strong agreement between experiment and traditional theory of 
the Casimir force was used to obtain more stringent constraints on the hypo- 
thetical Yukawa-type interactions predicted in high energy physics. Using the 
results of two experiments (this one and of [11]) the constraints on a as a func- 
tion of À were strengthened by up to a factor of 20 within a wide interaction 
range. 
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Table 1 

Casimir pressures and different contributions to them. See text for further discus- 
sion. 



Separation 


Impedance approach 




Plasma model 






z (nm) 


Pt 0) (P^ 


Pt~ 1] (Pa) 


Pi (Pa) 


p(i=0) /p Q ^ 
^ 3 (Fa) 


^3 (Pa) 




p(i>l) /p Q X 

(Pa) 


Ibi) 


-0.0715 


-1.0726 


-1.1441 


-0.0719 


-1.0430 


—1.1149 


-1.07715 


nan 

200 


-0.03834 


-0.47096 


-0.5093 


-0.03o4í 


-0.46335 


— 0.5018 


n /i ^70 /i 
— U.4/24 


nrri 

250 


-0.02046 


-0.20421 


-0.2247 


-0.02050 


-0.20260 


— U.zzol 


—0.2046 


nan 

300 


— U.U12ZU 


n i noi a 
— U.1U214 


— U.114o 


—U. 01222 


— U.lOloz 


n 1 i /i n 
— U.114U 


n 1 noQ 
—0.1023 


350 


-0.007858 


-0.056424 


-0.06428 


-0.007866 


-0.056425 


-0.06429 


-0.05648 


400 


-0.005358 


-0.033547 


-0.03890 


-0.005362 


-0.033616 


-0.03898 


-0.03357 


450 


-0.003817 


-0.021105 


-0.02492 


-0.003820 


-0.021178 


-0.02500 


-0.02111 


500 


-0.002816 


-0.013887 


-0.01670 


-0.002817 


-0.013948 


-0.01676 


-0.01389 


600 


-0.001659 


-0.006665 


-0.008324 


-0.001660 


-0.006703 


-0.008363 


-0.006667 


700 


-0.001059 


-0.003546 


-0.004605 


-0.001059 


-0.003569 


-0.004628 


-0.003547 



Table 2 

Casimir pressures and zero-frequency contributions to them in the alternative ap- 
proaches to the thermal Casimir force. 



Separation 
z (nm) 


Approach 1 [26,34,35] 
P[ l=0) (Pa) P 1 (Pa) 


Approach 2 [27] 
P 2 a=0) (Pa) P 2 (Pa) 


160 


-0.04836 


-1.1255 


-0.0967 


-1.1739 


200 


-0.0247 


-0.4971 


-0.0495 


-0.5219 


250 


-0.0127 


-0.2173 


-0.0254 


-0.2300 


300 


-0.00734 


-0.1096 


-0.0147 


-0.1170 


350 


-0.00462 


-0.06110 


-0.00924 


-0.06572 


400 


-0.00310 


-0.03667 


-0.00619 


-0.03976 


450 


-0.00217 


-0.02329 


-0.00434 


-0.02546 


500 


-0.00158 


-0.01547 


-0.00317 


-0.01706 


600 


-0.000917 


-0.007584 


-0.001835 


-0.008502 


700 


-0.000577 


-0.004124 


-0.001155 


-0.004702 
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Table 3 

Mean vàlues of differences between the theoretical and experimental Casimir pres- 
sures. See the text for the definition of Afi. 



z 


Atot, T nthooT" / \ noxnt / \1 

A [P iacoT {z) ~ P cxpl {z)\ 






A D 

APi 


A T~> 


(nmj 


(mFaj 


(mFaj 


(mFaj 


(mFaj 


(mFaj 


170 


17.2 


2.01 


13.0 


18.8 


-21.8 


180 


13.4 


-0.74 


7.54 


1 A A 

14.4 


-19.8 


200 


8.59 


-1.21 


5.3 


11.0 


-13.9 


250 


O O A 

3.34 


-0.31 


1.3 


7.09 


-5.66 


300 


1.59 


0.34 


0.6 


5.07 


-2.28 


350 


0.89 


0.38 


0.39 


3.58 


-1.05 


400 


0.63 


0.28 


0.20 


2.59 


-0.62 


500 


0.49 


0.11 


0.05 


1.37 


-0.22 


600 


0.46 


0.08 


0.04 


0.82 


-0.09 


700 


0.46 


0.02 


-0.01 


0.51 


-0.07 



Fig. 1. Schematic diagram showing the experimental setup. The different distances 
are defined in the text. 

Fig. 2. AFM images of the Au films deposited on the MTO and on the sphere. 
Topographies on the different panels are indicated. 10 x 10 /t/m 2 images of the film 
on the MTO before the experiment (a) and after it (b) are provided. Figures (c) 
and (d) provide similar information for the sphere, except that a 5 x 5 /im 2 area is 
shown. 

Fig. 3. Typical cross sections of the atòmic force microscope images of the Au 
coatings on the plate (a) and on the sphere (b). h denotes the height of the surface 
above the reference level defined in the text, and is plotted against the lateral 
position x. 

Fig. 4. Topographic heights hi on the sample as a function of the total sample area 
with heights h < hi + \ for the film deposited on the MTO (a) and the film deposited 
on the sphere (b). For both films the differences observed when analyzing different 
regions in the sample are negligible. 

Fig. 5. Dependence of the angular deviation 9 as a function of the applied voltage to 
the sphere. Data obtained at five different separations z between the metallic layers 
are shown. Data have been displaced vertically for the sake of clarity. 

Fig. 6. Electrostatic force F e i as a function of separation z for Va u — Vq = 0.22 V, 
0.25 V 0.35 V. For simplicity, the z-axis has already been corrected for D and 6. The 
lines are fits to the data using Eq. (3). 



43 



Fig. 7. Change in resonant frequency Af r = (u r — ü) r )/(27r) as a function of ampli- 
tude of motion of the MTO, obtained at a separation z = 300 nm. The non- linear 
behavior becomes dominant at large amplitudes, increasing the uncertainty in the 
frequency. Frequency shifts are measured with respect to the one obtained using the 
thermodynamic noise as a driving force, see text. Inset: Resonance curve obtained 
under the conditions indicated by the arrow. 

Fig. 8. Absolute value of the parallel plate Casimir pressure P cxpt as a function of 
separation z. The separation has been corrected taking into account D and 9. 

Fig. 9. The random A rand P cxpt (long-dashed lines), systematic A syst P expt 
(short-dashed lines) and total A tot P expt (sòlid line) absolute experimental er- 
rors of the Casimir pressure measurements versus separation within the intervals 
160 nm < z < 750 nm (a) and 160 nm < z < 750 nm (b). 

Fig. 10. The total relative error of the Casimir pressure measurements versus sep- 
aration. 

Fig. 11. The relative theoretical errors of the Casimir pressure computed in different 
approaches due to use of the proximity force theorem (short-dashed line) , variation 
of the tabulated optical data (long-dashed line), due to the experimental uncertain- 
ties of absolute separations (dotted line), and total theoretical relative error (sòlid 
line) versus separation. 

Fig. 12. The 95% confidence intervals (sòlid lines) and differences [i-4 heor — P^p*] 
(dots) versus separation within the intervals 160 nm < z < 250 nm (a) and 
250 nm < z < 750 nm (b) for one set of measurements computed in the impedance 
approach. 

Fig. 13. The 95% confidence intervals (sòlid lines) and differences [P4 heor — P^p*] 
(dots) versus separation within the intervals 160 nm < z < 250 nm (a) and 
250 nm < z < 750 nm (b) for all fourteen sets of measurements computed in the 
impedance approach. 

Fig. 14. The 95% confidence intervals (sòlid lines) and differences [P^ 2 Cor — P expt ] 
(dots) versus separation for theoretical approaches 1 [26,34,35] (a) and 2 [27] (b) 
incorporating different models of alternative thermal corrections. All fourteen sets 
of measurements as in Fig. 13 were used. 

Fig. 15. Constraints on the strength of the Yukawa interaction a versus interaction 
range À. Line la is obtained in this paper, line lb is from [11], line 2 follows from 
old Casimir force measurements between dielectrics [12]. Lines 3 and 4 are obtained 
from Casimir force measurements by using a torsion pendulum [2,13], and by means 
of an atòmic force microscope [6,16], respectively. The region of (a, À) plane above 
each line is excluded and below each line is allowed. 
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